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A  TRANSFER  MATRIX  FOR  MODULARLY- DRIVEN  ARRAYS 

By  H.  H.  Ding  and  J.  A.  Ward 


INTRODUCTION 

In  the  course  of  working  with  certain  sonar  transmitting  array  prob¬ 
lems,  we  must  determine  the  transfer  matrix  relating  element  head  veloci¬ 
ties  to  amplifier  input  signals.  We  require  this  matrix  to  predict  element 
head  velocities  as  a  function  of  a  given  distribution  of  amplifier  input 
signals  and  vice-versa  without  resorting  to  scale-model  or  prototype  array 
s  tud  Les . 


For  the  purposes  of  this  memo,  we  define  this  transfer  matrix  in  the 
following  manner  for  an  array  of  N  elements: 

[Si]  “  [Kij3  [vj]  (i’j  =  1>2>3 . 


where:  =  the  signal  into  the  amplifier  driving  element  r. 

v^  =  the  head  velocity  at  element  P. 

[K^  ]  =  the  N  x  N  transfer  matrix  ^ 


The  purpose  of  this  memo  is  to  present  the  development  of  this  trans¬ 
fer  matrix  for  several  types  of  modularly-driven  transducer  arrays.  The 
simplest  type  of  modularly-driven  array  is  that  wiih  only  one  element  per 
module.  The  development  of  ]  for  this  type  of  array  is  presented  first 
to  prepare  the  way  for  more  general  developments  for  several  types  of 
modularly-driven  arrays. 


We  would  like  to  emphasize  at  this  point  that  the  transfer  matrix 
used  in  this  memo  may  be  computed  utilizing  data  from  available  computer 


prog  i  unis.  No  significant  modification  of  existing  programs  is  required; 
rather  a  new  program  must  be  written  which  would  use  several  available 
programs  in  effecting  final  computation  of  the  [K^ }  transfer  matrix.. 
However,  for  large  arrays,  the  equations  developed  here  may  not  be  eco¬ 
nomically  feasible  without  further  simplification  and  development  of 
shortcuts  by  a  capable  computer  programmer. 

[K^]  FOR  THE  SIMPLE  CASE 

For  this  development,  we  will  assume  identical  elements  and  identical 
amplifiers  in  addition  to  assuming  one  element  per  amplifier.  For  the 
rth  element  of  an  N-element  array,  the  following  diagram  provides  the 
notation  used  here: 
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Figure  1:  The  rth  Module 

Knowledge  of  array  geometry,  medium  propagation  constants,  etc., 
allows  computation  of  mutual  coupling  impedances  for  some  array  shapes. 
For  our  N-element  array: 

[ Fi J  “  [Zij]  [v^J  (i,j  =  1.2,3, . . ,r . N) 

where:  F  =  force  at  the  head  of  the  rth  element, 
r 

v *  velocity  at  the  head  of  the  rth  element. 

[Zlj]  =  N  x  N  mutual  impedance  matrix, 
and,  obviously,  )  (F^)  =  )  <^zrJVj) 


'  1 


The  A-matrices  used  to  describe  the  amplifier  and  the  transducer 
e lenient  are  defined  as  follows: 
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Rewriting  the  transducer  equation  in  another  form: 
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By  means  of  the  fact  that: 
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Since  the  elements  are  identical,  we  may  write 


for  the  entire  array,  where: 


lziv 

] 

Izu] 

tvi 

1 

(i,j  -  1,2,3,. 

. . ,r, . . . 

where : 

did 
=  Zij  "  Zij 

(i  *  j) 

=  Z  +  Z 
ij  ec 

(i  -  j) 

2u 

=  Z'  +  Z, 
ij  ic 

(i  -  j) 

We  may  now  determine  [K^j ]  for  two  types  of  amplifier  drive. 


If:  S  =  El  ; 
r  in  ’ 


If:  S  =  r  ; 
r  in 


iKijj  * 

[Kij]  -  *1' 


(voltage  drive) 
(current  drive) 


For  voltage  driVe: 


[Si]  *  tEin]  *  a  [Ei]  +  8  [Ii] 
[S.3  =  (aA  [Z®  1  +  BC  IZ^])  [Vj] 

is,]  i  IK^J  (Vjl 


Thus:  [K8.]  =  aA  [Z^J  +  6C  [Z^J  (i,j  -  1,2,3 . 2-,...,N) 


Note  that  for  an  amplifier  that  behaves  as  a  perfect  voltage  source  (3  =  0)  , 
[K.Ej  ]  reduces  to  aA  (Z8^]. 

For  current  drive: 

tSi]  '  lIin]  “  y  lEi]  +  6  [li] 

[Si]  =  (yA  izj  1  +  6C  [Z^])  [Vj] 

>V  -  i^j1  (V 


♦ 
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(i,j  *  1,2,3 . r,..,N) 


[Kj  ]  =  YA  [zj  ]  4  6C  [z[  J 

Aga i i; ,  note  that  for  an  amplifier  that  behaves  as  a  perfect  current  source 
(y  -  0),  [K^J  reduces  to  6C  [Z^  ) 

E  I 

Once  [K_ ]  or  [K^  ]  has  been  determined,  it  is  a  straightforward  mat¬ 
ter  to  determine  [ S J  given  lvjJ  or  to  determine  [v  ]  given  [S^],  since 


and 


[Vj]  .  [SJ 

■si>  -  ("iji  (V 


(l.j  =  1,2,3,...,  N) 


Also,  inversion  of  [K^]  should  be  aided  by  the  symmetry  existent  fen  thic 
simple  case. 


Wo  can  expand  on  this  a  little  by  removing  the  assumption  that  all 
elements  are  identical  and  all  amplifiers  are  identical.  If  we  do  this, 
the  following  definitions  become  necessary,  for  the  rth  element  and  the 
t’th  amplifier: 


Finally;  [K^  ]  =  diag  (a1)  diag  (A^  ([Z^]  +  diag  (Z*c>  ) 

+  diag  (8i)  diag  (Cp  ( [Z±J ]  +  diag  (zjfi)  ) 

OC  [Kij]  =  dia§  (aiAi  +  ^i*  fZijJ  +  diag  (aiAiZec  +  6iCiZic) 
Similarly:  [k[  ]  =  diag  (y^  +  5^)  [Z^J  +  diag  (y^Z^  +  fi^Z^) 


This  form  readily  shows  that  the  off-diagonal  terms  of  [K^ ]  become 
unimportant  compared  to  the  on-diagonal  terms  when  either  Z  or  is 
made  large  enough.  One  further  simplification  can  be  shown  by  assuming 
the  amplifier  to  be  a  perfect  voltage  source  (or  current  source): 


Voltage  source:  8=0 


[K_J  =  diag  <o±Ai)  [Z±  J  +  diag  (a^) 


Current  source:  y‘=‘ 0 


[K^J  =  diag  (6^.)  [Z.^]  +  diag  (6^) 


i  E 

Thus,  if  Z  is  made  large  (B.  >>  A.),  the  diagonal  terms  of  [K . . ]  take 
ec  i  l  ij 

the  approximate  form: 


“i  (AiZit  +  V  '  “iAi  (Zli  +  Zec> 


where  Z  >>  Z.. 

ec  ij 

Similarly  for  ;  diagonal  terms  take  the  form: 
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and 


A  MORI-  GENERAL  CAS ^PARALLELED  ELEMENTS) 

an  N-eleineut  transmitter  array  consisting  of  G  modules, 
/on  by  a  single  amplifier.  The  transducer  elements  within  a 
idule  are  connected  electrically  in  parallel  and  are  Q  in  nunibe 
iry  from  module  to  module).  The  following  diagram  shows  the 
ao-l  elements  of  the  <yth  module  along  with  the  notation  used 
>•' o l  i  on 


Figure  2:  The  ;yth  Module  Couta  1  n;  i.g  Q  Elements 
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Note  that  the  number  of  elements  is  not  constrained,  nor  are  they 
necessarily  identical.  For  this  genera]  development,  we  will  not  assume 
identical  amplifiers  nor  will  we  assume  identical  elements. 


To  avoid  possible  snarls  in  notation  in  the  subsequent  development, 
we  choose  to  state  our  subscripting  convention  at  this  point. 


(i > j  =  1,2,3,. ...r, 
(k,m  1,2,3,. .. ,g, 
(n  =  1,2, 3, 

(p  =  1,2,3, ...  ,q, . . 


,  N) 

N  = 

total 

,  G) 

G  * 

total 

V 

^k  " 

total 

*»> 

total 

array  elements, 
array  modules, 
elements  in  kth  module, 
elements  in  mth  module. 


If  q  or  Q  appear  in  a  double  subscript  (e.g. :  gq  or  gQ)  it  will  be 

understood  that  q  or  Q  refer  to  elements  within  the  module  denoted  by  the 
other  subscript  (e.g.:  in  gq,  q  is  the  qth  element  of  the  £?th  module  and 
<?Q»  Q  i-S  understood  to  be  Q  ) . 


Thus,  we  can  now  write  [F^]  =  [Z^ ^ ]  [v^ ]  In  a  form  which  names  the 
modules  as  well  as  the  elements  and  which  is  more  suited  to  our  needs 
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here.  This  form  is  as  follows: 
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w !;>  •  r  ]/.  !  is  a  Q.,  x  Q  rectangular  matrix  in  general  and  where 

.!  j.  K  tn 

^'np',  is  an  N  x  N  square  matrix, 
km 

Moreover,  our  convention  can  be  extended  to  allow  representation  of 
i I  e  column  vector  components  and  single  square  matrix  elements  in  the 
i o 1  1  ov i ng  manner : 

F  is  the  force  at  the  q th  element  of  the  gth  module. 

;  7c/ 

Z  ,  is  the  mutual  impedance  relating  the  force  at  element  q  of 

module  g  to  the  velocity  at  element  s  of  module  h. 

v,  is  the  velocity  at  element  s  of  module  h. 

Hi-  J 


Note  that  the  ordering  of  subscripts  is  consistent  with: 
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(Q  =  Qg) 
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on  order) 


With  these  statements  in  mind,  we  write  (for  the  qth  element  of  the 
guti  module): 
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Because  the  elements  are  connected  in  parallel  within  a  single  module, 
we  write  (for  the  c/th  module): 


As  shown  previously,  the  mutual  and  self  impedances  are  written  as 
Eol lows : 


flF  1  1 

= 

riz  j  1 

J"(v  J  1 

L  n  kJ 

L  np  km  J 

L  P  mJ 

or,  for  short: 


[F.  ]  =  [Z.  ]  [v  ] 

kn  kmnp  mp 


and  this  again  leads  to:  (— — 1  (F  )  =  (-^— )  (E  Z  v  ) 

vgq  99  vgq  m.p  9*1 P  mp 


Rewriting  the  transducer  equation  for  the  q  th  element  of  the  cth 
module  gives: 
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==  A 

F  4-  B  v 

04 

04 

04  04  07 

I 

=  C 

F  +  D  v 

04 

04 

04  04  04 

E  =  A  (-  ^  ^  4  7pq )  v 

qq  qa  v  ec  gq 

'  gq 


I  =  c  (—<’  4  7?q )  v 

gq  gq  -gq  -  ,-v 


who n  : 


,  B 

7G‘l  d  _gq 

ec  A 

gq 

,  D 

7aq  d  _2 q 

ic  C 

gq 


K L  •  tl  iy:  E  =  A  (£  Z  v  +  Z^v  ) 

gq  gq  mjP  g^q p  “>p  ec  gq 


i  =  c  (r  z  v  +  zwy  ) 
gq  gq  mjP  «p  ®p  ic  <?q 


Wo  may  now  write  (for  the  entire  array  of  G  modules): 


[Ekn]  =  diae  (V  EZkmnp]  [Vmp] 

E 1 k n E  =  dlaS  (Gkn)  EZkmnp]  [Vmp^ 


whtre:  di  (A^)  ^  diag^  diag  (An) 

<C,p!  j  [l^W.] 


kmnp 

^krnnp 

Zkmnp 

e 

7. 

+ 

zkn 

kmnp 

kmnp 

ec 

z 

+ 

zkn 

kmnp 

kmnp 

1C 

Iff(k  ^  m  or  n  ^  p) 
Iff  (k  =  m  and  n  -  p) 
Iff  (k  =  m  and  n  =  p) 


Finally,  we  can  determine  the  transfer  matrix, 

both  the  instance  where  S  =  E,  and  where  S  =  I.  . 

in  in 


for 


[K, 


kmnp 


=  EKlnpE 
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[S^]  =  diag  (ai)  [Ei]  +  dlag  (3.^)  [Ii] 

or:  [S®]  =  diag  (a^)  [E^  +  diag  (3^  [Z Ikn] 

Since  we  wish  to  determine  the  relationship  between  the  head  veloci¬ 
ties  of  the  elements  and  the  signal  inputs  of  their  associated  amplifiers 
we  rewrite  the  above  equation  in  a  form  which  facilitates  the  desired 
development. 

lSaV  '  dl»*  «W  [tEn'k]  +  dl*E  <V  fEnVj 

where  [Se],  [E  ],  and  [II  ]  are  column  vectors  having  Q.  identical  com- 
n  n  n  n  K 

ponents . 


By  substituting  where  appropriate: 


[S 


knJ 


(diag  diag  (A^)  [<mnpl  +  diag  (6^)  diag  (C^) 

fez*  ]  1  \ 

L  a  °P  km J  ) 


x  [v  ] 

mp 


Thus : 


>:  ]  -  IKT  ]  [vm  ] 

kn  kmnp  mp 

or,  more  formally: 


ft se ]  1 
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[ke  ]  1 

rtv  1  1 

L  n  J 

l  np  km J 

L p  J 

whe  re : 


|jKnp^  j  =  dia8  ^iaB^an^  diag  (liag(An)^  |[Zn 
+  diag  (diag(Sr)^  diag  ^Jiag(Cn)^  ^ 
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np 
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The  development  in  this  case  is  nearly  identical  in  principle  to 
that  presented  in  the  previous  section.  The  results  will  be  presented 
here  without  comment. 


Q  [Snlk  = 


More  General  Modular  Arrays 

If  an  array  has  several  modules  wherein  some  elements  are  connected 
electrically  in  parallel  and  others  are  connected  in  series,  the  transfer 
matrix  for  the  array  may  be  developed  in  a  manner  similar  to  the  develop¬ 
ments  presented  in  this  memo.  We  do  not  include  such  a  development  here 
as  we  feel  that  it  would  be  of  very  limited  usefulness  at  the  present 
time . 
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